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Abstract 

For a hyperbolic Brownian motion on the Poincare half-plane H 2 , starting from a point of hy- 
perbolic coordinates z — (rj, a) inside a hyperbolic disc U of radius rj, we obtain the probability of 
hitting the boundary dU at the point (fj,a). For fj — > oo we derive the asymptotic Cauchy hitting 
distribution on cffl 2 and for small values of r\ and f\ we obtain the classical Euclidean Poisson kernel. 
The exit probabilities P z {T r)1 < T V2 } from a hyperbolic annulus in H 2 of radii r/i and 772 are derived 
and the transient behaviour of hyperbolic Brownian motion is considered. Similar probabilities are 
calculated also for a Brownian motion on the surface of the three dimensional sphere. 

For the hyperbolic half-space H n we obtain the Poisson kernel of a ball in terms of a series involving 
Gegenbauer polynomials and hypergeometric functions. For small domains in BP 1 we obtain the n- 
dimensional Euclidean Poisson kernel. The exit probabilities from an annulus are derived also in the 
n-dimensional case. 

Keywords: Hyperbolic spaces, Hyperbolic Brownian motion, Spherical Brownian motion, Poisson kernel, Dirich- 
let problem, Hypergeometric functions, Gegenbauer polynomials, Cauchy distribution, Hyperbolic and spherical 
Carnot formulas 

1 Introduction 

Hyperbolic Brownian motion has been studied over the years by several authors on the half-plane H 2 
and on the Poincare disc D 2 and more recently in the n-dimensional hyperbolic space (see, for example, 
Matsumoto and Yor [13], Gruet [TO], Byczkowski et al. [2] and Byczkowski and Malecki [3]). The 
hyperbolic half-space H™ is given by 

HP = {z = (x, y) : x e R" -1 , y > 0} 

with the distance formula 

u , / x 1 \\ z '- z \\ 2 

coshniz , z) — 1 H — - — . 

%yy 

The hyperbolic Brownian motion is a diffusion governed by the generator 
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(see, for example, Gruet QH]). Therefore the probability density p(xi, . . . , x n _\, y, t) of hyperbolic Brow- 
nian motion is solution to the Cauchy problem 



dp y 2 [ v-i d 2 p d 2 p \ (n — 2) dp 



Ot 2 \ ^-f dxf ' <9y 2 / 2 



subject to the initial condition 

n-l 

• • -j^n-lj^O) = ]J <5(j/ - 1). 



For our purposes it is important to express the generator ( |1.1[ ) in hyperbolic coordinates (77, a) = 
(77, ai, . . . , a„_i) as follows 



d 2 n — 1 d 

A n = — + ——— + — ^A Sn _ t 1.2 
dr\ z tanh 77 d?7 smh 77 



where Ag ii _ 1 is the Laplace operator on the (n — l)-dimensional unit sphere (see, for example, Helgason 
[TT] page 158 or Grigor'yan [8]). 

The aim of this paper is to study the hitting distribution on a hyperbolic sphere for a hyperbolic 
Brownian motion starting from an arbitrary point inside the sphere. Our work is related to the paper 
by Byczkowski et al. [2] where the Poisson kernel of half-spaces in i n , n > 2, is studied and the paper 
by Byczkowski and Malecki [3| where the Poisson kernel of a ball in the Poincare disc D", n > 2, is 
considered. 

The first part of our paper concerns the derivation of the Poisson kernel of a hyperbolic disc in H 2 by 
solving the Dirichlet problem 

tt~2 + <. v, #- + ■ u2 #r m m,a:n, a) = 0, < n < fj < 00, 

ar\ z tanh rj or} smh^ r] da z J W 1 *' / 1 ' ' ' (13) 

w(77, a; ?7, a) = 5(a — a), a, a € (— tt, 7r]. 

The interplay between Dirichlet problems and hitting probabilities in various contexts is outlined, for 
example, in Grigor'yan [5]. The explicit solution of (1.3) is 

1 cosh fj — cosh 

' ' ' 2ir cosh 77 cosh ^ — 1 — sinh?7sinh/7COs(a — a) 

and represents the hitting distribution on the hyperbolic circumference of radius fj for the hyperbolic 
Brownian motion starting at (j], a). 

We show that for fj — > oo the distribution (1.4) tends to the Cauchy distribution as was found by 
means of other arguments in Baldi et al. [Tj. 

The solution to the Dirichlet problem (1.3 1 is carried out by two different approaches. One is based 
on the direct solution of the hyperbolic Laplace equation and the second one is based on some integral 
representation of the associated Legendre polynomials. 

The derivation of the n-dimensional Poisson kernel for n > 2 is much more clumsy and the final 
expression is given as a series involving Gegenbauer polynomials and hypergeometric functions. A sub- 
stantial simplification in the calculations is obtained by applying a suitable rotation of the hyperbolic 
sphere so that the kernel can be expressed in terms of the hyperbolic distance 77 and the angle a% — a± 
between the geodesic lines with ends points (tj, a) and (fj, 6c). Its explicit form reads 

u(?7,Q!i;T7,ai) 
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= — — > ^ + 1 7 o^fCfc (cosfai - ai))sm™ ^(ai - aA 

n « fco^ n_2 / tanh fe 2f (fc,l-f;*; + i;tanh 2 2) fc V V " K h 

(1.5) 

where n > 2, < 7] < fj < oo, cti — fix € (0, t], = pfff) ^ s ^ ne sur f ace area of the n-dimensional 

Euclidean unit sphere, F(a, f3;j, x) is the hypergeometric function and C[ (x) are the Gegenbauer 
polynomials. 

The Poisson kernel of a ball in the hyperbolic disc O", n > 2, is obtained in Byczkowski and Malecki 
[3], formula (16) and must be compared with (1.5| above. 

Unfortunately formula (1.5| cannot be reduced to a fine form as (1.4). However, for sufficiently small 
domains, we extract from (1.5) the n-dimensional Euclidean Poisson kernel. 

Section 3 is devoted to the exit probabilities P z {T r)1 < T rj2 \ from a hyperbolic annulus of radii rji and 
7/2. We examine in detail both the planar and the higher dimensional case discussing also the transient 
behaviour of hyperbolic Brownian motion. 

In the last section the hitting probabilities on a spherical circle for a spherical Brownian motion 
starting from p — {$,(p) are considered. In particular the most interesting result here is that 

1 cos ?9 — cos ?9 — 

F p {B s (T,)edrt = -- = -dip, < < < tt, ^06(0,24 

Zir 1 — cos cos v — sin v sin v cos(</j — ip) 

2 Hitting distribution on a hyperbolic sphere in M n 
2.1 Two dimensional case 

We study here the Poisson kernel of the circle in the hyperbolic plane H 2 = {(x,y) : x E M.,y > 0} 
endowed with the Riemannian metric 

, 2 dx 2 + dy 2 
ds = , 

yZ 

and the distance formula 

2yy' 

We denote with 77 the hyperbolic distance from the origin O — (0, 1) of H 2 . The Laplace operator on H 2 
in cartesian coordinates reads 

(for a proof see, for example, Chavel [5 page 265). It is convenient to write the Laplace operator in 
hyperbolic coordinates (77, a) 

d 2 Id Id 2 

dr\ 2 tanh 77 drj sinh 2 77 da 2 

(for information on hyperbolic coordinates see Cammarota and Orsingher [4 ) . The relationship between 
hyperbolic coordinates (77, a) and the cartesian coordinates (x, y) is given by 

sinh rj cos a 

cosh rj— sinh rj sin a ' ^2 ^ 

y cosh r\ — sinh r\ sin a ' 



^~n2 \ tsnln <n pirn ni-nli^ ^ /")/n/2 ^ 



By exploiting (2.4), in the paper by Lao and Orsingher [12], the Laplace operator (2.3) is obtained from 
(2J2). 

We have now our first theorem. 
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Theorem 2.1. Let U = {(77, a) : 77 < 77} be a hyperbolic disc in H 2 with radius r\ and center in O, the 
solution to the Dirichlet problem 



1 d 2 



is given by 



dr\ 2 tanh r\ dr\ sinh 2 r/ da 2 

u(fj, a; fj, a) = 5(a — a), 



u(r),a;r),a) 



u(rj, a; rj, a) = 0, < 77 < 77 < 00, 
a, a € (— 7T, 7r], 

cosh 77 — cosh 77 



27r cosh 77 cosh 77 — 1 — sinh 77 sinh 77 cos (a — a) 



Proof 

Our proof is based on the classical method of separation of variables. We assume that 

u(?7, a; 77, a) = £'(77)9(0;) 
and we arrive at the following ordinary equations 

fe"(a)+/x 2 e(a)=0, 

1 sinh 2 77 E"(i]) + cosh?7 sinh 77 E'(n) — ^, 2 E(n) = 0, 
where /i 2 is an arbitrary constant. The first equation has general solution 

8(a) = Acos(pa) + Bsm(jia) 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



and becomes periodic with period 2ir for /i = m G N. The second equation necessitates some further 
treatment. We start with the change of variable w = cosh 77 which transforms the second equation of 
EM into 



(1 - w 2 )G"(w) - 2wG'(w) ^G(w) 

1 — w z 



The general solution to (2.101 can be conveniently written as 

m/2 



G(w) = Ci 



w + l 

w — 1 



Co 



w — 1 
w + l 



m/2 



= 0. 



m ^ 0, 



(2.10) 



(2.11) 



(see, for example, Polyanin and Zaitsev [Tl] Section 2.1.2, formula 233 for a = 1, b = —1, A = and 
fi = — m 2 ). From (2.11) we have that 



cosh 77 + I 
V cosh 77 — 1 



m/2 



( cosh 77 - 1 \ m/2 _ ( 
V cosh 77+I ) ~ Cl[ 



cosh 77 + I 
V sinh 77 



■a 



( cosh ?7 — l 
\ sinh 77 



(2.12) 



(2.5 1, so that we have 



We disregard the first term of (2.12| since our aim is to extract finite- valued and increasing solutions to 

(2.13) 



E( v ) = c 



( cosh 77 — 1 
\ sinh 77 



Ctanlr 



In light of (2.71, (2.9) and (2.13) we can write 



i(r), a; 77, a) = Q m (a)E m (n) = A + ^ [ A m cos(ma) + B m s\n{ma) 



m=0 



m— 1 



/ cosh 77 — 1 
\ sinh 77 



(2.14) 
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If we take the Fourier expansion of the Dirac delta function 

1 1 °° 

5(a — at) — - — I — / cos[m(a — all 



1 1 

2n 7r 



m— 1 

00 



[cos(ma) cos(ma) + sin(ma) sin(ma)], 



(2.15) 



by comparing (2.14) with (2.151 wc obtain the Fourier coefficients A m and B m so that we can write 
1 1 



i(n,a;r],a) = — + 
Zir 



2tt 
1 

2^ 



1 



lcos(ma) cos(ma) + sin(ma) sin(ma)l ( | [ 

V sinh?7 J V 

E, . _. . /cosh 77 — 1\ ~ m /cosh?? — 1 
cos(m(a - a)) . - _ — 
\ sinn 77 J \ smn 7/ 



cosh 77 — 1 
sinh 77 



00 r 



m— 1 



i(a-a) silm7 7 COsh^-l 



cosh 77 — ! sinh 



sinh 77 cosh 77 — 1 
cosh 77 — 1 sinh 77 



/ \ 2 

/ cosily — 1 
I sinh fj 



(cosh rj — 1 \ 
sinh r; J 



2tt 
1 



/ coshq-A 2 + / cosh^jV _ 2 cosM-l coshr,-! cos(a _ -x 
V sinn 17/ V sinn n y sinn r} smn r/ v ' 



(2.16) 



tanh 2 f - tanh 2 § 



27r tanh 2 2 + tanh 2 2-2 tanh § tanh \ cos(a - a) ' 
The expression in (2.16) can be substantially simplified by observing that: 

(cosh rj — l) 2 sinh 2 77 — (cosh 77 — l) 2 sinh 2 77 — (cosh fj — 1) (cosh 77 — 1)[2 cosh 77 — 2 cosh 77] 

and 

(cosh 77 — l) 2 sinh 2 77 + (cosh 77 — l) 2 sinh 2 77 — 2(cosh 77 — 1) (cosh rj — 1) sinh 77 sinh 77 cos(a — a) 
— (cosh 77 — 1) (cosh 77 — 1) [2 cosh 77 cosh 77 — 2 — 2sinh?7sinh77COs(a — a)]. 

In view of all these calculations we have that the hyperbolic Poisson kernel takes the form 

I cosh 77 — cosh 77 



u(n, a; 77, a) 



2-7T cosh 77 cosh 77 — 1 — sinh 77 sinh 77 cos (a — a) 



Remark 2.1. It is possible to obtain the expression ( 2.14 ) by means of an alternative approach as follows. 
We start from the associated Legendre equation 



(1 - z 2 )y"{z) - 2zy'(z) + 



1/(1/ + 1) 



1-z 2 



y(z) - 



(2.17) 



which coincides with (|2.10| for v — or v = — 1. In view of Gradshteyn and Ryzhik 9J formula 8.711.2, 

-dip, I argz| < 



the solution to (2.17) can be written as 

(_!)«• r(i/ + i) 



cos(mtp) 

7r T(h> — m + l)J [ Z + y/z 2 - 1 COSipY +1 
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If v = — 1 we have 



P™ x (cosh 77) = 



(-l) m 1 
7T r(-m) J 



cos(m(j))d(j) = 0, for raeZ. 



If 1/ = we have 



P m (cosh7y) 



1 



dcj) 



for m = 1,2, 



7r r(l - m) y cosh 77 + sinh 77 cos <j> 1^0, for 777 = 0,-1,-2, 



It follows that 



u(rj,a;fj,a) = [A m cos(ma) + B m sin(ma)]P ( 7 t (cosh?7) 



m,— — 00 




cos(m0) 



- ^ ^ r [An cos(ma) + B m sin (ma)] / 

7T ■* — ' 1 1 — 777 n 

m=-oo y 1 JtJ 

- V(-l) m L4 m cos(ma) +B m sin(ma)] / — . 

7T ^— ' ,/ cosh 77 + sinh 77 cos 



cosh 7/ + sinh 77 cos ■ 
cos(m^) 



(2.18) 



in the last step A n and B n include the multiplicative constant r( - 1 j_ m - 1 ■ Since we have 



>+*£('4^)"-<«*> 



1 



cosh 77 + sinh 77 cos < 



(2.19) 



by inserting (2.19) into (2.18) we get 

1/(77, a; 77, a) 



j 00 / 00 

- (-l) m [A m cos(ma) + B m sin(ma)] / cos(m^) 1 + 2 

™ — n ^0 V „ — 1 



1 — cosh 77 \ ™ 



J \ sinh 77 



cos(?70) d(f> 



2 00 00 / 1 cosh \ n f 1 
— ^2 (~ l) m [^m cos(ma) + 73 m sin(ma)] ^ ( . ^ - ) / cos(m</>) cos(?70)d0 



m=0 



= + [il m cos(ma) + f? TO sin(ma)] 



m— 1 



n=l 

/ cosh 77 — 1 
\ sinh 77 



sinh 7; 



and thus we retrieve (2.14). 



Remark 2.2. By applying the hyperbolic Carnot formula we note that it is possible to write the hyper- 



bolic Poisson kernel (2.6) in a new form. We construct a hyperbolic triangle with sides of length 77, 77 and 



77, and angle between the two sides of length 77 and 77 equal to 9 = a — a, see Figure [T] The hyperbolic 
Carnot formula 

cosh 77 = cosh 77 cosh 77 — sinh 77 sinh 77 cos(a — a), 



permits us to write (2.6) as 



u(r),a;r),a) 



1 cosh 77 — cosh 77 



2-7T cosh?) — 1 
where the dependence of u from a and a is hidden in 77. 



(2.20) 
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Figure 1: Hyperbolic triangle in H 2 with sides of length 77, 77 and 77. 



Remark 2.3. We observe that the hyperbolic Poisson kernel (2.6) is a proper probability law. In fact: 



• It is non-negative because, for 77 > 77, we have cosh 77 — cosh 77 > and by the hyperbolic Carnot 
formula 

cosh 77 cosh fj — 1 — sinh77sinh77cos(o! — a) — cosh 77 — 1 > 0. 



It integrates to one since it is well-known that 



/o a + b cos 9 yfa 2 - b 2 
where, in this case, a = cosh 77 cosh 77 — 1 and b = — sinh?7sinh?7. 



(2.21) 



Remark 2.4. The kernel appearing in formulas (2.6) and ( |2.20 ) represents the law of the position 



occupied by the hyperbolic Brownian motion {Sup (7;) : t > 0} on H 2 starting from z = (77,0) € H 2 when 
it hits for the first time the boundary dU of the hyperbolic disc U. In other words 

™ r„ , m % ,i 1 cosh 77 — cosh 77 

F z {B n , (T n € da} = — ; . - .'_ -da, ae[0,27r), 

Zir cosn 77 cosn 77 — 1 — sinn 77 smn 77 cos(a — a) 

where T n = inf{rj > : B M2 (t) € dU}, see Figure (2] 



Remark 2.5. For small values of 77 and 77 the hyperbolic Poisson kernel (2.6) is approximated by the 
Euclidean Poisson kernel 

1 l + f-q+j) 1 rf-77 2 

u(n, a; 77, a) ~ 5 = r - 

27r(l + ^)(l + ^)-l-, ? ^cos(a-a) 2tt ?f + rf - 2 77 fj cos(a - a) 

that represents the law of the position occupied by the Euclidean Brownian motion {B(t),t > 0} on K 2 
starting from a point z = (77, a) when it hits for the first time the boundary dU of the Euclidean disc 
U = {(77, a), 77 < 77} with Euclidean radius 77. This is a consequence of the fact that in sufficiently small 
domains of the Lobatchevskian space, the Euclidean geometry is in force. 
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Figure 2: Brownian motion on H 2 starting at {rj,a) Figure 3: Hyperbolic Brownian motion starting in- 
and hitting the boundary of the hyperbolic disc U. side the hyperbolic annulus A with radii 771 and 772. 



Remark 2.6. We also note that: 



• For 77 = formula (2.6) becomes the uniform distribution as expected. 



• For 77 — ^ 00 we have that 

71(77, a; a) :— lim 7/(77, a; fj, a) 



1 1 

27r cosh 77 — sinh 77 cos(a — a) 



(2.22) 



In view of (2.19), the limiting distribution (2.22) can also be written as 

/ cosh T] — 1 



u(rj, a; a) 



1 



1 



1 



27T cosh 77 — sinh 77 cos(a — a) 2tt 



00 

n=l 



\ sinh 77 



cosn(a — a) 



We note that u represents the hitting distribution of the hyperbolic Brownian motion, starting at z = 
(7/, a), on the horizontal axis 9H 2 = {(f), a) : fj = 00} = {(x,y) : y = 0}, see Figure [i] We observe that 
the boundary <9H 2 represents the point at infinity of H 2 . We can write the 'hitting' probability on <9H 2 
in the following form 



^-BbpCT,,,,) e da} = - 1 — ^ — : ^da. 

Zir cosn 77 — sinn 77 cos a cos a — sinn 77 sin a sin a 



We write now the distribution (2.23) in cartesian coordinates. In view of (2.4) we have that 

-2 ly 2_ 1 



X 

y 



sinh 77 cos a, 



tan a = 



x" 



2x 



(2.23) 



(2.24) 



The first relation is an immediate consequence of (2.4) and for a proof of the second equality, see Cam- 

(2.25) 



marota and Orsingher [5]. From (2.1) and (2.24) it follows that 

tana 



sinh 77 sin a = y cosh 77 — 1 



x 2 + y 2 - 1 



Vl + tan 2 1 



2y 
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Letting 77 — > 00 we note, in view of (2.4), that for a point (x,y) £ dM 2 it holds that 




Formula (2.26) implies that x — cos a = x\f\ — cos 2 a and this leads to the following relations 

7,2 



2.7: 



1 



cos a 



1 



1 



(2.26) 



(2.27) 



In view of (2.24), (2.25) and (2.27) and since da = 1 _^- 2 dx, we can write u(r],a;a)da in cartesian 
coordinates as tollows 



u(x, y; x)dx 



1 



1 



2^ x 2 +y 2 + l 

1 



x 2x 
y 1+.T 2 



2y 



-11-s 2 1 

1+x 2 



:dx 



2y 



-dx 



7T (x 2 + y 2 + 1)(1 + x 2 ) - Axx - (x 2 + y 2 - 1)(1 - x 2 ) 

1 y 

7r x 2 (a; 2 + y 2 ) — 2xx + 1 

1 y 

7T 



d.c 



x^Jx 2 + y 2 - , x 

V 1 +v 



-dx 



x z -\-y 



1 



x 2 -\-y 2 



x 2 +y 2 



7 da;. 



(2.28) 



Formula (2.28) says that the probability that the hyperbolic Brownian motion starting at (x,y) £ H 2 

hits the boundary of H 2 at (x, 0) is Cauchy distributed with scale parameter y' = x 2+ y 2 and position 

parameter x' — 2 °L 1 depending on the starting point. In particular, if the hyperbolic Brownian motion 
x -\-y 

starts at the origin O of H 2 , we obtain a standard Cauchy. We note that (2.22) can be viewed as a Cauchy 
density in hyperbolic coordinates. 



Remark 2.7. In view of formula (2.28), we also note that the probability that the hyperbolic Brownian 
motion starting at z = (x, y) — (77, a) £ H 2 hits <9H 2 at (x, 0) is equal to the probability that a Euclidean 
Brownian motion starting at z' — (x 1 , y') — (77, a') hits the cc-axis at (x, 0), where z and z' have the same 
hyperbolic distance rj from the origin but a' = —a, see Figure [5] In fact 



cosh rj 



tana' = 



(x 2 +y 2 ) 2 



y 

(x 2 -\-y 2 ) 2 



y 



1 



2y 



x 2 +y 2 



(x 2 +y 2 ) 2 



y 

(x 2 +y 2 Y- 



- 1 



2/y 

1 — x 2 — y 2 



cosh?7, 



2x 
x 2 -\-y 2 



2.T 



tana. 



Formula ( 2.28 1 is in accordance with formula (1.2) in Baldi et al. p]. In this paper the hitting distribution 
on the horizontal axis, for the hyperbolic Brownian with horizontal and vertical drift, is obtained from 
the hitting distribution on the horizontal lines H a = {(x, y) £ H 2 : y — a > 0} when a — >• 0. 



Remark 2.8. The Poisson kernel (2.6) can be conveniently written also in cartesian coordinates by 
exploiting the relations (2.24), (2.25) and the hyperbolic distance formula 

x 2 

cosh rj = 



2y 



1 
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We have that 

u{x,y;x,y) — — x2+y2+1 x2+y2+1 ^ x x x -2 +y 2_ 1 x 2 +y -2_ 1 

~ 2y 2y y y 2y 2y 

_ 1 {x 2 +y 2 + l)y - {x 2 + y 2 + l)y 
7r x 2 + y 2 + x 2 + y 2 — 2yy — 2x5 

= 1 (x 2 + y 2 )y - (a; 2 + y 2 )y + y - y 
7T (x - x) 2 + (y - y) 2 

In the special case where y = the previous expression becomes 

1 (l + x 2 )y 
u(x,y;x,0) = i — j 

and thus multiplying by 1 _^- 2 we get the Cauchy density as expected. 




Figure 4: Brownian motion on H 2 starting at (t),o) Figure 5: Hyperbolic Brownian motion starting at 
and hitting the boundary of the hyperbolic plane. z and Euclidean Brownian motion starting at z' . 



2.2 Multidimensional case 

Let H™ = {z = (x,y) : x £ E" _1 ,y > 0} be the n-dimensional hyperbolic plane, n > 2, with origin 
O = (0, . . . , 0, 1) endowed with the Riemannian metric 

2 da^ + .-. + diz^+dy 2 
as = - 



\z' -z\\ 2 



and the distance formula 

coshnfz', z) = 1 + 

with r) :— i](0, z). The Laplace operator on H n in cartesian coordinates is given by 

Q2 q2 \ Q 
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(for a proof see, for example, Chavel [3] page 265). The Laplacian A„ in hyperbolic coordinates (77, a) = 
(r),ai, . . . ,a„_i), reads 

d 2 n — 1 d 1 .„ 

drf tanh 77 d?7 S mh 77 

where Ag ii _ 1 is the Laplace operator on the (n — l)-dimensional unit sphere (see, for example, Helgason 
HO page 158). 

In Lem ma|2.1| we evaluate the hyperbolic Laplacian of the distance 77 in H™. This result permits us, 
in Theorem |2.2| to determine the hyperbolic Laplacian of a smooth function f(r]). The statement of this 
result is given, for example, in Davies [7] page 117 without proof where coshp must be replaced by cothp. 

Lemma 2.1. For z = (x,y) and 2' = (x' ,y') in H™ we have that the hyperbolic distance 77(2,2') is a 
solution of 

A ( ') n — 1 

tanh 77(2, 2') 

Proof ,22/2 

Since coth(arcosh(x)) = ^ r ^,_ l and 77(2, 2') = arcosh W x ~ x [L±g +v | we to prove that 
A„77(2, z ) = in — 1) — . 

V[lk-^ll 2 + (2/ + y') 2 ] [lk-^ll 2 + (2/-y') 2 ] 



In fact we have 



#^,*') = Hx-^ + t/ 2 -, 2 = , (2.30) 

dy Vy /[ \\x - x'\\ 2 + {y + y 1 ) 2 ] [ \\x - x'\\ 2 + {y - y') 2 } 

d 2 . I Ice — .t'II 2 + y' 2 — y 2 

; r)(z,z) 



dy 2 ' v 2 \J\ W - x '\\ 2 + {y + y') 2 }[\\x- x'\\ 2 + {y~ y') 2 } 

4\\x - x'\\ 2 [\\x ~ x'\\ 2 + y 12 + y 2 } 

+ ^[\\x-x>\\ 2 + {y + yi) 2 Y - A\ 2 + {y -WY' 

On the other side, for i = 1, . . . , n — 1, we have 

9 / / \ _ 2{Xj — Xj) , . 

~ y/[\\x- x'\\ 2 + (y + y') 2 } [ \\x - x'\\ 2 + {y - y') 2 } ' 



d 2 , 2 

^r (z ' z) ~ 7TF ^'ll 2 + (y + y') 2 ] [ ¥ - x '\\ 2 + (2/ - y'H 



A(xi - x'i) 2 [\\x - x'\\ 2 + y 2 + y' 2 



y/[\\x- x'\\ 2 + (y + y') 2 } 3 [\\x- x'\\ 2 + {y - y') 

v- 1 , A 2(n-l) 

S^i y/\\\x- x'\\ 2 + (y + y') 2 } [ \\x - x'\\ 2 + (y - y') 2 } 
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k\\x - x'|| 2 [ ||x - x'|| 2 + y 2 + y' 2 



y/[\\x- x'\\ 2 + (y + y') 2 ] 3 [ \\x - x'\\ 2 + (y - y') 2 ] 3 



So, finally, we obtain that 

A n r)(z,z') = y 2 



" _1 <9 2 d 2 
.i=i 1 



d 

(71-2)77—77(2,2') 
dy 



11 



2(n - l)y 2 

y/[ \\x-x'\\ 2 + (y + y') 2 } [ ||x-x'|| 2 + {y-y 1 ) 2 } 

\\x — x'\\ 2 + y' 2 — y 2 
+ ^[\\x-x'\\ 2 + {y + y>) 2 ] [\\x-x'\\ 2 + {y-y'f 
| {n - 2 )[\\x-x'\\ 2 + y' 2 -y 2 } 

^[\\x-x'\\ 2 + {y + y') 2 \ [\\x-x>\\ 2 + {y-y>) 2 

(n-1)- 



\x — x'\\ 2 + y' 2 + y 2 



y/[\\x- x'\\ 2 + {y + y') 2 } [ \\x - x'W 2 + (y - y') 2 



In view of Lemma |2.1| the following theorem holds: 
Theorem 2.2. If f is a smooth function on R, it holds that 

n- 1 



Proof 
We have 



A n /(r?) = f'(v) + 



A n f(r,(z,z')) 



= y 



d 2 



y 

drj 2 ' 
drj 2 



,i=i 

n-l 

E 



d 2 f ( dr) 



drj 2 \dx 



E 1 

i=l 
n-l 

E 

.*=! 



dr) 
dxi 

dn 

dxi 



d 2 f df . 
orf or) 



tanh.77 



/'(*?)• 



(n-2)y—f(r)(z,z')) 
dy 



5/ AV &l (dvV + dl 
dr) dx 2 J dr) 2 \dy J dr) dy 2 



dr) 

dy 

dr) 
dy 



df 
df) 

d£ 
dr) 



' 1 d 2 r) d 2 r) 



E 



dxf dy 2 



- (n-2)y 
-(n-2)y 



df drj 
dr) dy 

dr) 



dy 



since, in view of formula (2.30) and (2.31), it holds that 

n-l 



E 



dr) 

dx,. 



,2 12 



4| |ac - x'\\ 2 y 2 + [ \\x - x'\\ 2 + y' 1 - y 
y 2 [ | \x — x'\ | 2 + (y + y') 2 ] [ | \x — x'\ \ 2 + (y — y') 2 ] y 2 



From Lemmal2.ll we obtain the final result 



dr) 2 dr) dr) 2 tanhrj dr) 



We denote with {B^n(t),t > 0} the hyperbolic Brownian motion on H™ with starting point z = 
(r),a) € EP where a = (aj., . . . , oe n -i) G [0) 7r ] n_2 x [0, and we assume that z is inside the n- 
dimensional hyperbolic ball U with hyperbolic radius fj (see Figure |6|. 
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We are interested in obtaining the law of the position occupied by the hyperbolic Brownian motion 
on HP when it hits the boundary dU for the first time. 

Since the Laplace operator is invariant under rotations (see, for example, Helgason [11] Proposition 
2.4), without loss of generality we can assume that the starting point is z = (77, oc\, 0, . . . , 0) and the 
process hits the boundary of the ball U at some point z — (fj, cl\, 0, . . . , 0), where ot\ — 6t\ is the angle 
between the vectors z and z. For a function on the (n— l)-dimensional unit sphere SVi-i depending only 
on one angle 9 we have 



8 



1™- 2 de 



*1 

de 



if 



n-2 8 



d0 2 tan 6 89' 



(2.32) 



In view of (2.291 and (2.32 1 we have that the hitting distribution on dU is obtained from the solution of 
the following Dirichlet problem 

n-l d 1 1 / d 2 
nh 2 rj \ da f 



_il ( 

dr] 2 ~ r tanh rj dr) ^~ sinh 2 rj V da 2 ~ r tan ai dai 

u(fj,ax;fj,ai) = 6(ai - ai), 



u(r), ot\\r), ot\) = 0, < 77 < 77 < 00, 
ai — ai & (0, 7r]. 



(2.33) 



Theorem 2.3. The solution to the Dirichlet problem {2.33) is given by 
u(i7,ai; 77,61) 



bin 



E 



2k 



, 2 2^ 



\ tanh fc | f ( fc; i„ t;fc+ | ;tanh ^ j c( ^ = 



-. , , . , fcosfai — ai)) sin™ 2 (ai — ai), 

fc ^ o Vn-2 J tanh fc 2 F (fc ; i_| ;fc+ | ;tanh 2 2) fe v v x " V1 ly ' 



(2.34) 



where n>2,0<i]<rj<oo and a.\ — a,\ € (0, n] . 
Proof 

As in Theorem |2 . 1 1 our proof is based on the method of separation of variables. We assume that 

u(77,Q!i;?7,ai) = Q(ax)E(r)). 

E(rj) 



Since we have that 



B(a 1 )E"(r 1 ) + e(a 1 )^-J-E'(r,) 



tanh 77 



sinh 2 r\ 



tanai 



0, 



there exists a constant /u 2 such that 



e"(ai) + (n - 2) cot a x &{a x ) + // 2 9(ai) = 0, 

sinh 2 77 E"{rj) + (n — 1) cosh 77 sinh 77 E'(r)) — ^ 2 E(-q) = 0. 



(2.35) 



The first equation in (2.351 can be reduced to the Gegenbauer equation. With the change of variable 
uj = cos ai and for /i 2 = k(k + n — 2), we obtain 



(1 - oj 2 )G"{lu) - (n - l)wG'(w) + k(k + n- 2)G{w) = 0. 



(2.36) 



The Gegenbauer polynomials 2 (w) satisfy (2.36) (see, for example, Polyanin and Zaitsev [2] S. 2. 11-4 
or Helgason [11] page 16) and this implies that 

9(ai) = AC ( k ^ } (cos ai). (2.37) 
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We transform the second equation of (2.35) 



sinh 2 T)E"(ri) + (n - 1) cosh rj sinh r]E'(r]) - pL 2 E(r]) = 



(2.38) 



into a hypergeometric equation. The first step is based on the change of variable £ = tanh | . We have 
that 

d 1 d 



dr; 2cosh 2 fdC' 
_cP_ 1 d 2 

dr/ 2 



sinh 3 d 



4 cosh 4 I dC 2 2 cosh 3 f dC 



By taking into account that sinh 77 = 2 sinh 3 cosh 5 and that cosh r/ = 2 cosh 2 § — 1, equation (2.38) 
becomes 



4 sinh - cosh 

2 2 



_JL d^ 

4cosh 4 2dC2 2 cosh 3 ? dC 



sinh ^ d 



E(0 



+ (n- 1)2 sinh ^ cosh 'j- ( 2 cosh^ - 1 



2 ?7 



2 cosh 2 2 dC 



E(C) - ti 2 E(0 = 0. 



And since 



2 sinh 



3 n 



sinh 3 



+ (n - 1) aUU1 1 (2cosh 2 - - l) = tanh- -2sinh 2 - + (n - 1) (2cosh 2 - - 1 



cosh ^ 



' cosh ^ 



2 



2 ?7 



?7 

tanh — 
2 

77 

tanh — 



1 + (n - 2) (2 cosh 2 I - 1 

1 , , 1 + tanh 2 § 

1 + (n-2) 2^ 

1 — tanh 5 



we can write equation (2.38) as 



tanh 2 -E"(0 +tanh- 
2 2 



1 + (n - 2) 



1 + tanh 2 § 
1 - tanh 2 a 



^'(0-^(0 = 0, 



that is 



C^"(C) + c 



1 + (n - 2) 



1 + C 2 
1-C 2 



B'(C) - M 2 S(C) = 0. 



We now assume that 



(2.39) 



E(C) = C fe /(C 2 )- 



Since 



^(CHfcC^W^C^V'tC 2 ), 

E"(0 = fc(fc - i)C fe - 2 /(C 2 ) + 2fcC fe /'(C 2 ) + 2(fc + i)C fe /'(C 2 ) + 4C fe+2 /"(C 2 ) 
- fc(fc - i)C fc " 2 /(C 2 ) + 2(2fc + i)C fc /'(C 2 ) + < k+2 f" (C 2 ), 

by replacing (|2.40| and (|2.4l[) into (|2.39|) (with pi 2 = k(k + n - 2)) we have that 



(2.40) 
(2.41) 



*(* - i)C fe /(C 2 ) + 2(2fc + i)C fc+2 /'(C 2 ) + 4C fc+4 /"(C 2 ) 



14 



1 + (n - 2) 



i + C 2 



i-C 2 

and with obvious simplifications we have that 

k(k - 1)/(C 2 ) + 2(2fc + i)C 2 /'(C 2 ) + 4C 4 /"(C 2 ) 
i + C 2 ~ 



[k( k f(( 2 ) + 2C fc+2 /'(C 2 )] -k(k + n- 2)C fc /(C 2 ) = 



1 + (n - 2) 



i-C 2 



[fc/(C 2 ) + 2C 2 /'(C 2 )] ~k(k + n-2)f(( 2 ) 



4C 4 /"(C 2 ) + 2C 2 



2(k + 1) + (n - 2) 



i + C 2 



i-C 2 



/'(C 2 ) + 2(n-2)fc r ^/(C 2 )=0. 



After some additional manipulations we arrive at the following equation 

C 2 (i - C 2 )/" (C 2 ) + [k + \ - (k + 2 - \ ) C 2 ] /' (C 2 ) + * (I - i) /(C 2 ) = o. 



(2.42) 



Equation (2.42| coincides with the hypergeometric equation 

i(l - t)/"(t) + [7 - (a + + l)t]f'(t) - a(3f(t) = 

for t = ( 2 , a = k, /3 = 1 — § and 7 = fc + |. In view of the position F(£) = C fc /(C 2 ) an d C — tanh 2, we 
conclude that a solution to (2.38) is given by 



F(??) = tanh fc I F (k, 1 



n n 2 V 

— ; k H ; tanh — 

2' 2' 2 



(2.43) 



Equations (2.37) and (2.43) imply that 



4(77, ai; 77, o-i) = ^ E k {rj)@ k {ai) 

00 

= ^A fc tanh fe | f(m 



fc=0 



-;£;+-; tanh 2 



2 ?7\ 



C) 2 ; (cosai). 



In order to determine the coefficients A k by applying the boundary conditions we have that 



u(?7,ai; 77,0:1) = 5(ai - ai) = ^4 fc tanh fc - F [k, 1 

fc=0 



-;fc+^;tanh 2 |)cf" ) (cosa 1 ). 



By multiplying both members by C m 2 (cosai) sin™ - a.\ and then integrating we have that 

/ 71-2 

E(ol\ — ai)C m 2 (cos ai) sin™ - oti dai 



= A k tanh fc I F (fc, 1 - J fc + tanh 2 |) / 2 2) (cos o<i)ci" 2 ^(cos a x ) sin" -2 ot x dai 



fc=0 



= A m tanh 7 ™ - F m, 1 ;m + -; tanh" 



because the functions C^(x) form an orthogonal system on the interval x € (—1,1) (see Gradshteyn 
and Ryzhik [5] formula 7.313). This implies that 



2 77^ 7r2 3- ™r(m + n - 2) 

n !( m + ^)r(=^) s 



(cos ai) sin™ 2 cci 



m!(m+^)r(^) ; 



tanh m f F(m,l- f;m+f;tanh 2 f ) 7r2 3 -™r(m + n - 2) 
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We finally obtain 



u(r},ati;r),ati) 



-^"sin"" 2 ^ 



— — V , ,; 1 \i~ . i " r~ — ' ^ cr* - ' (cos ax)cr s ~' (cos a x ) 

2 3 -"tt ^ T(fc + n - 2) tanh fe 2 F (jfe, 1 - f ; fc + f ; tanh 2 §) fe V U 



jfc! (jfc 



rt-2 
2 



tanh fc § F (jfc, 1 - f; jfc + f; tanh 2 £ 



TO/ 



By rotational invariance and since C fc 2 (1) 

u(rj,ai;fj,ai) 

r(^) 2 sin"" 2 (a 1 -ai) 

fc=0 



(™ + fe 3 ) ' * ne ^ as * expression reduces to 



2 3 -™(n-3)!7r 



E * 



-2\ tanh fc 2 F(k,l 



n . u , n 
2 ' 2 



2 / tanh fc § F(fc,l-f;fc+f;tanh 2 f) 



;tanh 2 f n^-. 



We arrive at formula (2.341 by observing that 



2 r2„_i 

2 3 -"(n- 3)!tt ~ n-2 fl n 



where Q n — r(n/2) ^ s suriace area of the n-dimensional Euclidean unit sphere. In fact, since 
2 2 - n y^T(n - 1) = r(f )r(f - §), we have 



'r^-2^ 



1 



(f-l) : 



2 3 -"(n-3)!7r 2 3 "™(n- 3)!tt (n _ ^ 2 

2 r(|) r(f) 



n-2 



1 



r(f) s 



2 3 -"(n - 3)!tt _ iy 



2 r i 



i 



n-2 Vtt 2 2 -"V^r(n-l) n-2 r(f-i) 

2 fin-l 



n — 2 f2„ 

and this concludes the proof of the theorem. 



Remark 2.9. We note that for small values of n and n we obtain the Euclidean Poisson kernel. In fact, 
since tanh \ ~ §, Cj n) (i) = 2nt, C^i) = 1 and kC { k n) (t) = 2n[tC i k n _ + 1 1) (*) - C^ 1} (t)], we have that 



E 



2A: 



\ tanh fe 2F(fc ; i_| ; A; + | ; t a nh 2 2 
1 ' fcg " 



J \n-2 J tanh fc | F(fc,l - §;jfc 



..^If^^CcosCar-ar)) 

2 ] idilll 2 J 



E 

2 



n-2 



2fc 

n-2 4 

OO 

E* 



A =2 



n-2 



E 

fc=0 



(cos(ai — ai)) 



C fc 2 (cos(ai — ai)) + (n — 2)- cos(ai — ai) 



Cf 2 ^cosfoi-ai)) 



cos(ai - ai) ^ ( - J C fc ^ 1 (cos(ai — ai)) 

fc=2 
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E 



k=2 
1 

V 



1 / V 

C^ 2 (cos(ai - ax)) + r cos(ai - ai) + ? 



cos(ai — ai) 



2 oo 

E 

fe=0 



c£ 5) (cos(ai -ai)) 



C7< S) (cos(«i - ai)) + C^\ C os{ax - ax)) 



k=0 



cos(ai - ai) + 



2 1 



2?7 
277 



cos(ai — ax) 



4) 


2 


- COS(Ct!i 






.V 






__2 \ 






V 2 ) 



- cos(ai - ax) - 



f} 2 



ill 



2ry , _ , 7? 2 
cos(ai - ai) + ^ 



i] 



2r, 



i - \ V 
cos(ai - ax) + ^ 



2 \ 1 



1- % 



1 - ^ cos(ai - ax) 



Remark 2.10. The kernel ( 2.34 1 represents the marginal, with respect to a.2, ■ . ■ , a n -x, of the distribution 
of the position occupied by the hyperbolic Brownian motion {B^n (t),t > 0} starting from z = (r),a) G H n 
when it hits for the first time the boundary dU of the n-dimensional hyperbolic hypersphere of radius fj. 
For z = (77, 0), such distribution is given by 



,{Bn^m) £ da} = J2 



fc=0 



tanh fe § F (k,l - f + f ;tanh 2 \ 



2k 



n - 2 ' J tanh fe § F (jfe, 1 



(2.44) 



where n > 2, 77 < 77 , ai G [0, tt) is the angle between z and z, and 

1 



/(a) = -J- sin" 2 ax sin" 3 0:2 ... sin a n -2 



is the uniform density on S n -x- 



Remark 2.11. We observe that (2.44) is a proper probability law. In fact: 

• The non negativity is due to the non negativity of solutions of Dirichlet problems with non-negative 
boundary conditions. 

• It integrates to one, in fact 



/ ... / / ¥ z {B mn (T n ) eda} 
Jo Jo Jo 

\ tanh fc 2F(fc,l-f;/c+f;tanh 2 f) /•» 
/ tanh fc f F(fc,l-f;A:+f;tanh 2 |) J 



s 1, 



E 

fe=0 



2k 
n-2 



f 71-2 \ 

C k 2 (cosai) sin" _2 aidai 



fi«-i F (0, 1 



2 > 2 



tanh 



2 T7\ /»7T 



ft n 



F(0,1- f;f;tanh 2 f) 7 



C, 



(^) 



(cosai) sin" 2 ax dax 
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since, if k > 0, we have 



f C[ n) (cos 9) sin™" 2 9 dO = 
Jo 

(see Gradshteyn and Ryzhik 9J formula 7.311.1) and F(0,/3;j;z) = 1, Cq(x) = 1, 

1 n-V 



sin™" 2 9 d6 = B 







Remark 2.12. We also note that: 



• For 77 — >■ (i.e. when the starting point is the center of the hyperbolic hypersphere) formula (2.441 
becomes the uniform distribution on SVi-i as expected. 

• For 77 —t 00, since tanh | — » 1 and F(a, /3; 7; 1) = p^^^"!^] if 7 > a + /3 (see Gradshteyn and 
Ryzhik [9 formula 9.122.1), we have that 

lim V z {B^(T n ) £ da} 



fc=0 



Remark 2.13. Byczkowski et al. in [2] provide an integral formula for the hyperbolic Poisson kernel of 
the half-space H a = {(x, y) £ H™ : y > a} for n > 2, a > 0, and show that for a — >■ 0, it converges to the 
Cauchy-type distribution 

r(n - 1) / j/ 



Tr^r(^) vr + F 

3 Exit probabilities from a hyperbolic annulus in HP 

3.1 Two dimensional case 

Suppose the hyperbolic Brownian motion {B$p(t),t > 0} starts at z = (r/,a) £ H 2 inside the hyperbolic 
annulus A with radii < 771 < 772 < 00 

A = {(??,«) : 771 < 77 < 772} 

(see Figure [3| . We define the hitting times 

T m = inf{* > : V (0,B B 2(t)) = 77J, i = 1,2, 

and T = T m AT, 2 . In the next theorem we evaluate the exit probabilities ¥ z {T m < T^}. Since these are 
given in terms of harmonic functions on the annulus A, they are closely related to the Dirichlet problem. 

Theorem 3.1. Let {B$p(t) : t > 0} be a hyperbolic Brownian motion starting at z — (r),a) £ A. The 
following result holds true 



log tanh % — log tanh % 
{ ^ <r - } = logtanhf-logtanhf * < * < ^ 
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Proof 

Since the probability in (3.1) is spherically symmetric we are lead to study the solution v : (771 , 772) 
to the Laplace equation involving only the radial part: 



d 2 



1 d 



dr] 2 tanh rj di] 



v{rj) = 



subjected to the boundary conditions v{r}x) = 1 and ^(772) = 0. With the change of variable w — cosh 77 
we immediately get 

(1 - w 2 )K"(w) - 2wK'(w) = 0, (3.2) 

whose general solution is 



K(w) = C\+C 2 log 



w — 1 



w+ 1 



(see, for example, Polyanin and Zaitsev [H] . Section 2.1.2 Formula 233 for a = 1, b = —1, A = and 
\i = 0). It follows that 



ufa) = Ci + C 2 log f C ° S ^ ■ ] ) =C X + C 2 log tanh J. 

\ cosh 77 + 1/ 2 

By imposing the boundary conditions we get 

11(772) - v(rj) log tanh ^ - log tanh 2 



(3.3) 



^(772) — log tanh ^ — log tanh ^ ' 



Starting from (3.1) and letting 772 go to infinity we have that Theorem 3.1 leads to the following 
corollary. 



Corollary 3.1. For any z — (77, a) outside the hyperbolic disc of radius 771 and center O, we have 

cosh 77+1 j log tanh I 



V z {T m < cx)} = 



log / cosl "i-i 



log 



(cosh 771 —1 
cosh 771+I 



log tanh % 



771 < 77. 



(3.4) 



It is possible to show with simple computations that the functions in (3.1) and (3.4) are genuine 
probabilities since they vary in (0, 1). 



3.2 Multidimensional case 

It is possible to generalize the exit probabilities from a hyperbolic annulus to the case of the n-th dimen- 
sional hyperbolic Brownian motion. 

In order to evaluate the exit probabilities from the hyperbolic annulus A in H™ , with hyperbolic radii 
771 and 772 with 771 < 772, we are interested in obtaining a solution v n : (771, 772) — > K to the radial part of the 
hyperbolic Laplace equation in H™ . We have proved in Lemma |2.1| and Theorem |2.2| that it is equivalent 
to solve 

r d 2 71 - 1 d 
d?7 2 tanh 77 d?7 



(3.5) 



In what follows we will assume that 



c(n,0) = 1, 



c(tt,, k) = 



(n - 3){n - 5) ■ ■ ■ (n - 2k - 1) 
(7i-2)(n-4)---(n-2fc-2) : 



k = 1, 



n — 3 
2 
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Theorem 3.2. For a hyperbolic Brownian motion {B^(t) : t > 0} started atz = (n, a) € A, we have that 
For n = 3, 5, 7, . . . 



eS(- 




1 c(n, fc) 


cosh r]2 


cosh 77 


+ (-!)" 


(n-3)!! 


log 


tanh ^y- 


sinh 7l_2fc_2 772 


s | n j 1 Ti-2fc-2 ^ 


(n-2)!! 


tanh S 


n-4 

Sfc=0 (~ 


-l) k - 


_1 c(n, Ar) 


cosh r/2 


cosh 7/1 


+ (-1) 


(n-3)!! 


lo£ 


tanh ^ 


s'mh n ~ 2k ~ 2 r/2 


sinh 71-2 ^ -2 771 


(n-2)!! 


' tanh ^ 



^ z\J-ni J- n 2 

For n = 4, 6, 8, . . . 
Proof 

The general solution to equation (3.5) is given by 

Unfa) = Ci +C 2 

For n = 2m + 1, m = 1,2,... we have 
1 



n — 3 
Efc=0 (" 


-l) k ~ 


1 c(n, k) 


cosh r/2 


cosh 7/ 


s inh n— 2fe— 2 r/2 


sinh 71-2 ^ -2 77 


71-3 

Sfc=0 (~ 


l) k ~ 


1 c(rx, A;) 


cosh r/2 


cosh r/i 


sinh 71-2 ^ -2 772 


sinh n— 2fe— 2 771 



(3.6) 



1 



z^—drj. 



sinh 77 



fnfa) =C\+C 2 



C!+C 2 



sinh 2m 7? d?? 



cosh 1] 
2m- 1 



1 y^, llfc -i 2 fc (m-l)(m-2)--.(m-fc) 1 

m l " -1 (2m - 3) (2m - 5) ■ ■ • (2m - 2k - 1) sinh 2 " 1-2 * -1 77 



sinh 2 ™- 1 7? 



smn ?y 



(3.7) 



fc=0 



(see Gradshteyn and Ryzhik [9] formula 2.416.2). 

For n = 2m + 2, m = 1, 2, . . . we have 

1 



u n (»7) = Ci + C 2 



Ci + C 2 



sinh 2 "^ 1 / 77 



cosh 77 
2m 



1 



sinh 2m 77 



EM) 



k=l 



fc _ 1 (2m-l)(2m-3)---(2m-2fc + l) 1 

2 k (m- l)(m - 2) ■ ■ ■ (m - k) sinh 2 '" -2 ' 8 1 



Ci + C 2 



fc=0 



sinh r ' 



(n-2)!! 



(3-8) 



(see Gradshteyn and Ryzhik ,9| formula 2.416.3). With computations analogous to those performed in 
the two dimensional case, we obtain the statement. I 



From this it follows immediately that: 
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Corollary 3.2. For z = (rj,a) outside the hyperbolic ball in H™ with radius r\\ and center in O, we have 
that 



For n = 3, 5, 7, . . . 



\{T. m < <x} 



For n = 4, 6, 8, . . . 



2 {T r)1 < cx)} 



Ego(-l) fc c(n,fc) sinh c °fe ni + (-I)? ^Ef 



(n- 


-4)!! 


(n- 


-3)!! 


(n- 


-4)!! 


(n- 


3)!! 


(n- 


2)!! 


(n- 


3)!! 



cosh rj 

sinh 7/ 



cosh 7)1 
sinh 7]i 



log tanh I 



t/i < r). 



log tanh ^ 



m < V- 



Remark 3.1. For the space H 3 formula (3.6) takes the simple form 

7/1 <Tj < 7/2, 



coth 7/2 — coth 77 



coth 772 — coth 7/1 ' 



and for 7/2 — > 00 yields 



1 — coth 7/ 
j{T r)1 < 00} = < 1. 



COth 7/! 

This shows that there is a positive probability that the hyperbolic Brownian motion never hits the ball 
of radius 7/1 . 



Remark 3.2. We note that for small values of 77 we have ^"^p 1 ^ ~ ^ and log tanh | ~ log 77. From (3.3) 
([3~7| and (pDH) it follows that 



sinh p r/ rj p 

d + C 2 log 77, ifn = 2, 

C1 + C277 2 -", if n = 3,4, 5. 



This means that, for sufficiently small domains, we obtain the exit probabilities of Euclidean Brownian 
motion from an annulus: 



P z {T m <T n2 } 



log r/ 2 -log 7)1 ' ' 



v 2 



— , if 77 = 3, 4, 5 . 



(3.9) 



Remark 3.3. It is important to note that for a planar hyperbolic Bownian motion the probability that 
the process goes to infinity before hitting the hyperbolic circle of radius rji is strictly less then one 



z {T Vl < 00} = 



log tanh 3 
log tanh \ 



< 1 



while it is well known, see (3.9), that for a planar Euclidean Brownian motion it holds that 

F z {T rn < 00} = 1. 

Hyperbolic Brownian motion is, in fact, transient for every dimension n > 2 as stated in Grigor'yan [5] 
Proposition 3.2. 
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4 Hitting distribution on a hyperbolic circle in D 2 

It is possible to obtain analogous results by considering a different model of the hyperbolic plane. In 
particular in this section we consider the Poincare disc model D 2 instead of the half-plane model H 2 . 
The half-plane H 2 can be mapped onto the disc D 2 = {(r,0) : r £ [0,1), 9 £ (— it, tt]} by means of the 
conformal mapping / : H 2 — ► D 2 such that 



/(*) 



iz + 1 



(4.1) 



The x-axis of H 2 is mapped onto <9B 2 while the origin O = (0, 1) of H 2 is mapped into the origin O = (0, 0) 
of D 2 . An arbitrary point z = (x, y) £ H 2 is mapped into a point Q = (r, 9) £ D 2 such that 



2r cos 9 



y 



1+r 2 — 2r sin t 

l-r 2 
1+r 2 — 2r sin 6 



(for details see Lao and Orsingher [H]). In view of (2.4 1 and (4.2) we have 

2r cos 9 



y 



— sinh r\ cos a = 



l-r 2 



Since we have that 



cos a = cos f 



and 



sinh?7 = 



2r 
l-r 2 ' 



for 9, a £ (— 7T,7r], we easily arrive at 



_ coshrj-l _ / coshrj-1 _ V 
sinh t; y cosh 77+! 2' 



(4.2) 



(4.3) 



The hyperbolic metric and the distance formula in D 2 become 



ds 2 = dr 2 

(l-r 2 ) 2 ' 



d(0, Q) = log 



1 + r 
l-r' 



By means of (4.2) the hyperbolic Laplacian in (2.2) is converted into 

2 2 \ld_ f d_\ ]_£?_' 

and the Dirichlet problem for the hyperbolic disc U = {(r, 9) : r < f } in D 2 reads 



(1 - r 2 ) 5 



r dr 



j_jr_ 

r 2 QQ2 



u(r,0;f,0) = 8(9-9), 



u(r, 9; r, 9) = 0, < r < r < 1, 
9,9 £ (— tt.tt]. 



(4.4) 



Since (1 — r 2 ) 2 > 0, we can derive the Poisson kernel related to the Dirichlet problem (4.4) from the 
Euclidean case: 



u(r, 9; r, 9) 



2ir r 2 + r 2 — 2rr cos(6* — i 



(4.5) 



coordinates. In fact, in view of (2.16) and (4.3), formula (4.5) immediately follows 



Alternatively it is possible to obtain formula (4.5) from the Poisson kernel in H with a change of 
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The Poisson kernel in (4.5) represents the law of the position occupied by the hyperbolic Brownian 



motion {B B 2 (t) : t> 0} on starting from Q = ( 
dll. We have 

¥ Q {B B 2{Tr) e d6} = ' 



when it hits for the first time the boundary 



-de. 



2-7T r 2 + r 2 — 2rr cos(0 — 



We note that for r — > 1 we have 



u(r, 9; 9) := lim u(r, 9; r, 9) = 

r->-i 2ir 1 



1-r 2 
r 2 — 2r cos(9 — 



r < 1. 



(4.6) 



Again u(r, 9; 9) represents the law of the position occupied by the hyperbolic Brownian motion in B 2 
when it hits for the first time the boundary of the hyperbolic disc U with hyperbolic radius that goes to 
infinity. Result (4.6 1 is stated in Helgason [TT] page 34. For the n-dimensional case see Byczkowski and 
Malecki [3] formula (16). 



In view of (4.3), we can write (4.6) in hyperbolic coordinates as follows 

2 

" 

n=l 

1 



1 



1 



2tt 1 + r 2 - 2rcos(6» - 9) 



1 
1 



1 OO 



r cosn f — t) = 



1 
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1 ^> / cosh 77 — 1 
7r ^— ' V sinh n 

n—l v ' 



cos n( 



2ir cosh 77 — sinh r\ cos(9 



which coincides with (2.22). On the other side it is well-known that under the conformal mapping (4.1) 
the Poisson kernel (4.6) takes the form of the Cauchy distribution as it is shown in formula (2.28). 

Since (1 — r 2 ) 2 > 0, the exit probabilities from the hyperbolic annulus A — {(r, 9) : rj < r < r 2 } are 
easily derived from the Euclidean case. If the hyperbolic Brownian motion starts at Q = (r, 9) G A, we 
have 

log r 2 - log r 



\{T ri <T r2 } 



Letting r-i — >• 1 in (4.7) we obtain that 

P Q {T ri < co} 



log r 2 - log T\ 
logr 



logri 



< 1, 



< n < r < r 2 < 1. 



< n < r < 1. 



(4.7) 



5 Brownian motion on the surface of a three-dimensional sphere 

The surface S of the unit-radius three dimensional sphere is a model of the elliptic geometry if geodesic 
lines are represented by great circles. We specify the position of an arbitrary point p 6 S with the couple 
($, ip) of spherical coordinates where •§ £ [0, n] and (p E [0, 2tt). 

If U = {(#, ip) : 1) > is the surface of a spherical cap on S with center in the south pole, the 
Dirichlet problem on the surface of the sphere S reads: 



d 2 , 1 

dt3 2 taniS 3i3 



^+ 1 



sin 2 ^ c)^ 2 



0, < i9 < < TT, 
€ [0,2tt). 



Assuming that u(i), tp; d, ip) = T(t?)i p (yj) we immediately arrive at the following ordinary equations 

\F"{ V )+p 2 F{p) = Q, 



sin 2 ?? T"(i?) + cost? sintf T'(^) - p? T(i?) = 0, 



(5.1) 
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with fi € M. With the change of variable w = cos#, in the second equation of (5.1), we arrive at equation 
(2.10) with general solution (2.11). Therefore, for /i = m e N, the general solution to the second equation 
of (5.1 ) can be written as 



= Ci 




(5.2) 



We restrict ourselves to the increasing component of (|5.2|) so that we have 



m=0 



[Am cos(mip) + B m sm(rmp)] 




By imposing the boundary condition u(i?, <p\ ip) = 5(<p — <p) and in view of (2.15), we finally obtain that 

1 , , , , / / 1 — cos $ / 1 + cos ■d 



u{-9, ip; i?, if) 



1 

2n 

1 



m— 1 



cos(m(c/2 — 93)) 




1 — cos i9 



1 



1 — COS 1? 1-fCOS 1? 

1+cos $ 1— cos d 



2w 



1 



1 



1 — cos $ 1+cos $ 
1+cos $ 1— cos (9 



1 — cos ^ 1+cos d 
1+cos i9 1 — cos $ 



COs(</? — <y9) 



cos •& — COS 1? 



2vr 1 



cos i9 cos $ — sin $ sin $ cos(ip — </?) 
Remark 5.1. Since for spherical triangles the following Carnot formula holds 

cos d = cos i? cos 1? + sin $ sin § cos(<p — ip) , 



(5.3) 



we can rewrite (5.3) as follows 



1 cos ■& — cos -d 
2tt l - cos d 



(5.4) 



Remark 5.2. We note that the Poisson kernel (5.3) is a proper probability law. In fact 



In view of (5.4) and observing that d < d we have that (5.3) is positive. 



Applying ( 2.21 ) with a = 1 — cos ■& cos ■& and b = — sin ■§ sin d we obtain that ( 5.3 ) integrates to one. 



For t? = 0we obtain from (5.3) the uniform law, while for $ = § we get 

7T 1 COS "0 

B( ^ ; 2^ ) = Sl- S mtfcc^-^)' 

Remark 5.3. Let {Bg(t) : t > 0} be a Brownian motion on the surface of the three dimensional sphere 
S with starting point p = (i9, tp) G S* (see Figure[7|. The kernel in (5.3) represents the law of the position 
occupied by the spherical Brownian motion when it hits for the first time the boundary of the spherical 
cap U. 

In order to obtain the exit probabilities of {B$(t) : t > 0} from a spherical annulus A = {(d,(p) : 
$2 < $ < $1} with center in the south pole of S, we consider the solution v : ($2, $1) — >• M to the Laplace 
equation involving only the radial part 



d-d 2 



1 d 

tan i? dd 



= 0. 
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Figure 6: Brownian motion on H 3 starting at (rj,a) Figure 7: Spherical Brownian motion starting at 
and hitting the boundary of the hyperbolic ball (0, ip) and hitting the boundary of the spherical disc 



With the change of variable w = cos0 we arrive at equation (3.2). With calculations analog to those 
performed in the proof of Theorem |3.1| we get 



In particular for 02 



log 


COS #2—1 


- log 


COS $— 1 


COS $2 + 1 


cos $+1 


log 


COS $2 — 1 


- log 


cos $1 — 1 


COS $2 + 1 


COS $1+1 



^ formula (5.5) reads 



02 < ■& < 01- 



(5.5) 



? p {T tfl < TV} 



log 


COS $ — 1 


COS $ + 1 


log 


COS $1 — 1 


COS $1+1 



log I tan I J 
log I tan ^r* 



< < 01. 



Remark 5.4. We note that replacing formally with id it is possible to extract from (5.3) and (5.5) the 
Poisson kernel and the exit probabilities obtained for the hyperbolic plane, namely (2.6) and (3.1). This 
is because H 2 can be viewed formally as a sphere with imaginary radius. For small values of we obtain 
instead results analogous to these obtained for the Euclidean Brownian motion. In fact, for sufficiently 
small domains, Euclidean geometry is in force. 



References 

[1] Baldi, P., Casadio Tarabusi, E., Figa-Talamanca, A., Yor, M.: Non-symmetric hitting distributions 
on the hyperbolic half-plane and subordinated perpetuities. Rev. Mat. Iberoam., 17, 587-605 (2001) 

[2] Byczkowski, T., Graczyk, P., Stos, A.: Poisson kernels of half-spaces in real hyperbolic spaces. Rev. 
Mat. Iberoam., 23, no. 1, 85-126 (2007) 

[3] Byczkowski, T., Malecki, J.: Poisson kernel and Green function of the ball in the real hyperbolic 
spaces. Potential Anal, 27, 1-16 (2007) 



25 



[4] Cammarota, V., Orsingher, E.: Travelling randomly on the Poincare half-plane with a Pythagorean 
compass. J. Stat. Phys., 130, 455-482 (2008) 

[5] Cammarota, V., Orsingher, E.: Cascades of Particles Moving at Finite Velocity in Hyperbolic Spaces. 
J. Stat. Phys., 133, 1137-1159 (2008) 

[6] Chavel, I.: Eigenvalues in Riemannian Geometry. Academic Press, New York, London (1984) 

[7] Davies, E., B.: Heat kernels and spectral theory. Cambridge Tracts in Mathematics, 92. Cambridge 
University Press, Cambridge (1990) 

[8] Grigor'yan, A.: Analytic and geometric background of recurrence and non-explosions of the Brownian 
motion on Riemannian manifolds. Bull. American Math. Soc. 36, 135-249 (1999) 

[9] Gradshteyn, I., S., Ryzhik, I., M.: Table of Integrals, Series, and Products. Fourth edition. Academic 
Press, New York-London- Toronto (1981) 

[10] Gruct, J-C: A note on hyperbolic von Miscs distributions. Bernoulli, 6, no. 6, 1007-1020 (2000) 

[11] Helgason, S.: Groups and Geometric Analysis. Academic Press, New York, London (1984) 

[12] Lao, L., Orsingher, E.: Hyperbolic and fractional hyperbolic Brownian motion. Stochastics, 79, no. 
6, 505-522 (2007) 

[13] Matsumoto, H., Yor, M.: Exponential functionals of Brownian motion. I. Probability laws at fixed 
time. Probab. Surv. 2, 312-347 (2005) 

[14] Polyanin, A., D., Zaitsev, V., F.: Handbook of Exact Solutions for Ordinary Differential Equations. 
Second edition. Chapman & Hall/CRC, Boca Raton (2003) 



26 



